Landau level spectra and the quantum Hall effect of multilayer graphene 
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The Landau level spectra and the quantum Hall effect of ABA-stacked multilayer graphenes are 
studied in the effective mass approximation. The low-energy effective mass Hamiltonian may be par- 
tially diagonalized into an approximate block-diagonal form, with each diagonal block contributing 
parabolic bands except, in a multilayer with an odd number of layers, for an additional block de- 
scribing Dirac-like bands with a linear dispersion. We fully include the band parameters and, taking 
into account the symmetry of the lattice, we analyze their affect on the block-diagonal Hamiltonian. 
Next-nearest layer couplings are shown to be particularly important in determining the low-energy 
spectrum and the phase diagram of the quantum Hall conductivity, by causing energy shifts, level 
anti-crossings, and valley splitting of the low-lying Landau levels. 

PACS numbers: 73.22.Pr 81.05.ue,73.43.Cd. 



I. INTRODUCTION 

Since the isolation of graphene flakes [l[ , the chiral na- 
ture of quasiparticles in monolayer and bilayer graphene 
has been observed in a range of phenomena including 
the integer quantum Hall effect [2|-|^] , Klein tunne ling [5| - 
H|, weak localization [§-[3, and photoemissio n Il7l-l20l . 
Generally speaking, the effective mass models [2rH30l j of 
monolayer and bilayer graphene have been very success- 
ful in describing these phenomena. The low-energy band 
structure of bila yer graphene is composed of a pair of 
parabolic bands [28-3Q), and is distinct from monolayer 
which has a Dirac-like linear dispersion [22|, [Hj]. In a 
ma gne tic field, the Landau level structures of monolayer 
[23H27I ] and bilayer [28|, [3l[ differ in the degeneracy at 
zero-energy, and quantum Hall plateaus appear at differ- 
ent filling factors accordingly [2|-|4|. 

In graphene multilayers with three or more layers, N > 
3 where N is the number of layers, the effective mass 
model is much more complicated than in monolayer or 
bilayer graphene f28,:29], or, even, bulk graphite [32| . On 
the one hand, there are more relevant parameters than 
in monolayer or bilayer because of the presence of next- 
nearest layer couplings while, on the other hand, a lack of 
translational invariance in the direction of layer stacking 
means that the number of basis states in the model is 2N, 
not four as in the Slonczewski-Weiss-McClure model of 
bulk graphite [32j . Although multilayers with a moderate 
number of layers, N > 10, arc thought to be similar to 
bulk graphite [33|, the properties of few- layer graphene 
(typically N>3) are distinct [H [H-|50|. 

Here, we consider ABA-st&cked (Bernal) multilayer 
graphene, and analyze its Landau level spectrum and 
quantum Hall conductivity in magnetic fields. In the 
lattice structure in Fig. Qla) and (b), the layers, each 
having two inequivalent atomic sites, A and B, on a hon- 
eycomb lattice, are stacked so that half of the atomic sites 
have a counterpart directly above or below in the adja- 
cent layers, referred to as dimer sites, whereas half the 
sites do not have such a partner (non-dimer sites). The 
effective mass model is characterized by intralayer cou- 



pling 70, nearest interlayer couplings 71, 73 and 74, and 
next-nearest layer couplings 72 and 75 |2ll. l22l. I2H [23. [32j . 
By comparison with experiments, it has been possible to 
determine the values of relevant parameters including in- 
tralayer coupling [^, [3|, [5ll 1521] and, in bilayers, interlayer 
couplings (l7l. [53H58j . 

While it is possible to numerically diagonalize an ef- 
fective Hamiltonian with 2N components ;'>.''> .'>3- this 
doesn't always shed light on the roles of the different 
parameters. Nevertheless, it was noticed that the multi- 
layer bands form groups of bilayer-like parabolic bands 
near the corners of the Brillouin zone called valleys [59[ 
with, in odd-iV layers, an additional pai r of monolayer- 
like bands with linear dispersion (29l . [33j . Subsequently, 
it was realized that a partial diagonalization of the ef- 
fective Hamiltonian could be performed in order to write 
it in block diagonal form [36H38] , as illustrated in Fig. []] 
(c) . Each block on the diagonal describes four bilayerlike 
bands with a given effective mass (labeled 'b' or 'B') or, 
in odd-iV layers, two monolayerlike bands (labeled 'M'). 

The decomposition into block diagonal form is based 
upon the eigenstates of a one-dimensional tight-binding 
chain in the stacking direction (perpendicular to the lay- 
ers) with nearest-neighbor hopping [29]. If the next- 
nearest layer couplings, 72 and 75, are neglected, the 
decomposition is exact: the diagonal blocks (and their 
corresponding bands) are completely separate from each 
other. The matrix elements associated with the next- 
nearest layer couplings appear within each block of the 
decomposed Hamiltonian, but they also couple separate 
blocks [as shown by the W blocks in Fig. [2(c)] and tend to 
hybridize their bands [381 ] . In this paper we calculate the 
Landau levels with the band parameters fully included, 
and show with the aid of decomposition that the next- 
nearest neighbor couplings generally account for the en- 
ergy shifts, the level anti-crossings, and valley splitting of 
the low- lying Landau levels, which significantly influence 
the phase diagram of the quantum Hall conductivity. We 
focus on trilayer, fourlayer and fivelayer, Fig. [TJ as rep- 
resentative examples. 
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FIG. 1: (a) Atomic structure of ABA-stacked multilayer 
graphene. (b) Side view of trilayer (left), fourlayer (center) 
and fivelayer (right) lattices, showing the mirror planes for 
odd- TV layers and the inversion center for even- TV" layers. Hor- 
izontal solid lines indicate 70 intralayer coupling between A 
(white) and B (black) atoms, vertical solid lines represent 
71. (c) Schematic form of the partially diagonalized Hamilto- 
nians with M, b, B labeling monolayerlike, light bilayerlike, 
and heavy bilayerlike blocks, W represents coupling between 
them arising from next-nearest layer couplings, 72 and 75. 



II. THE EFFECTIVE-MASS MODEL OF 
MULTILAYER GRAPHENE 

A. The effective-mass Hamiltonian 

To describe the electronic properties of Bernal-stacked 
multilayer graphene, we use an effective-mass model 
with the Slonczewski-Weiss-McClure parameterization of 
graphite [32| ■ The low energy spectrum is given by states 
in the vicinity of the K% point in the Brillouin zone, where 
£ = ±1 is the valley index. If \Aj) and \Bj) arc Bloch 
functions at the Kg point, corresponding to the A and B 
sublattices of layer j, respectively, then a suitable basis 
is \A 2 ),\B 2 ); \A N ),\B N ). In this basis, 

the Hamiltonian of multilayer graphene with N layers 
[H H HI H, HI in the vicinity of the K(. valley is 
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Here, the in-plane momentum operator is n = £p x + ip y , 
and p = {p x ,Py) — — iftV + eA with vector poten- 
tial A. The diagonal blocks, Eq. describe nearest- 
neighbor intralayer coupling, and V, Eq. describes 
nearest-neighbor layer coupling, where 71 is the inter- 
layer coupling between dimer sites. Parameter A' repre- 
sents the energy difference between dimer sites and non- 
dimer sites, and thus it only exists for N > 2. It is 
related to the band parameters as A' = A — 72 + 75 . The 
Fermi velocity of monolayer graphene is v = -\/3a7o/2?i, 
and other velocities are defined as v$ = ^/ia^z/2h and 
V4 = v3a74/2ft. Matrix W, Eq. (|4]), describes coupling 
between next-nearest neighboring layers, and it only ex- 
ists for N > 3. Parameters 72 and 75 couple a pair of 
non-dimer sites and a pair of dimer sites, respectively. 



B. Decomposition to an approximate block 
diagonal form 

The decomposition of the Hamiltonian, Eq. ([T]), into 
an approximate block diagonal form uses a unitary trans- 
formation based on the eigenstates of a linear chain of 
atoms in the z-direction 136 381: 
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where 



7T raw 
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(8) 



1/2 (m = 0) 
1/V2 (m^O). 

Here j — 1,2, ••■ ,N is the layer index, and m is the 
block index which ranges as 



m 



1,3,5,- • • ,N - 1, iV = even, 
0,2,4,- •• ,N- 1, N = odd. 



(9) 



Obviously f m (j) is zero on even j layers, while Qmjj) is 
zero on odd j layers. The basis is constructed 36|, |38( by 
assigning f m (j), g m (j) to each site as 
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where X = A or B. A superscript such as (A, odd) 
indicates that the wave function has a non-zero amplitude 
only on \Aj) sites with odd j's. 

In order to write the Hamiltonian Eq. (JT]) in terms of 
the basis states Eq. (fTU|) . we group the basis of block to 

as u m = {|0l A '° dd) >, |0l B '° dd) ), |^° vcn) ), |4 B ^° vcn) )}. 
Then, the block matrix between different m's may be 
written as 

Hm'm = U^j-ffUm — 7~L(\m)5m'm 

),(n) 

with 
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An, = 2 COS K, 
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which, barring the diagonal terms, is equivalent to the 
Hamiltonian of monolayer graphene. 

We stress the role of the symmetry of the lattice and 
note that the even-odd effect, with respect to the num- 
ber of layers, goes further than the absence or presence of 
monolayerlike bands in the band structure. The lattice 
of odd- A multilayers obe ys mirror reflection symmetry 
(x, y, z) — > (x,y 7 —z) [H, [4j| [6(| (mirror planes 

for trilayer and fivelayer graphene are shown in Fig.[TIb)), 
and thus the eigenstates can be classified by parity with 
respect to the reflection: the parity of the wavefunction 

N — TTl — 1 

of the subband to is given by (— 1) 2 , so the group of 
to = 2, 6, 10, . . . and that of to = 0, 4, 8, . . . have opposite 
parities. Since eigenstates with different parity cannot be 
mixed by terms in the Hamiltonian that preserve lattice 
symmetry, off-diagonal blocks connecting diagonal blocks 
with different parity are identically zero even in the pres- 
ence of next-nearest layer couplings [H, [Hf . We actu- 
ally see that the coupling matrix W mm < between blocks 
having different parities indeed vanishes, as illustrated in 
Fig. [TJc) for trilayer and fivelayer. 

Even-A multilayers lack mirror reflection symmetry, 
however, so that 72 and 75 mix every diagonal block, 
as shown in Fig. [TJc) for fourlayer graphene. Instead, 
the lattice of even-A multilayers obeys spatial inversion 
symmetry (x, y, z) -> (-x, —y, -z) jH, HE S3, HE 
(an inversion center for fourlayer graphene is shown in 
Fig. [U center). Unlike mirror reflection, inversion sym- 
metry transforms electronic states between valleys, and, 
even in the presence of a magnetic field, this ensures de- 
generacy of the electronic spectra at different valleys [45[ . 
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The diagonal matrix Hmm is equivalent to the Hamil- 
tonian of bilayer graphene [28j with nearest-layer cou- 
pling parameters multiplied by A [IE |33|, I36I4381 ] , and on- 
site asymmetric potential described by W mm . The off- 
diagonal block, W for to ^ to', has not been explicitly 
obtained before: it appears only when coupling between 
the next-nearest neighboring layers, 72 and 75, is non- 
zero. The block W is diagonal, where 72 only connects 
pairs of non-dimer sites, and 75 only connects pairs of 
dimer sites, of these effective bilayer-like blocks. 

The case of m = is special in that g m (j) is identically 
zero, so that only two basis states {\4> A ,odd - 1 ), |</>Q B ' odd ' 1 )} 
survive in Eq. (fTU|). The matrix elements associated 
with the two missing basis states should be neglected 
in Eqs. (|12I13I) . Specifically, the matrix for the to = 
block written in the two component basis is [38j 
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C. Reduced low-energy Hamiltonian 



As we show below, mixing between blocks is particu- 
larly important in the vicinity of level crossings. Level 
crossings apart, a good approximation to the spectra over 
a broad range of energy may be obtained by neglecting 
the off-diagonal blocks [Hj]. Taken alone, the bilayerlike 
block Hmm for to ^ describes four bands |28|], two 
split off by energy ±A 7 i at the point and two near 
zero energy. The split bands can be viewed as a bond- 
ing and anti-bonding pair created by the relatively strong 
interlayer coupling 71 between dimer sites (A, even) and 
(B, odd). For low energy, e -C A71, it is possible to de- 
rive a reduced Hamiltonian for the bilayerlike block de- 
scribing an effective hopping between non-dimer sites (A, 
odd) and (B, even) by using a Schrieffer- Wolff transfor- 
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This reduced Hamiltonian is approximately valid at low 
energy {e,vp} < A71. 

When W3 and W4 are neglected, the eigenvalues of Ti. mm 
at zero magnetic field are 
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with A = A r , 



a 



This gives nearly- 



parabolic conduction and valence bands centered at en- 
ergy (a + /3)72/2 with an energy gap | (a — /?)72 1 between 
them. The gap always vanishes for even-layered graphene 
since a mm — /3 mm holds for all m when N is even. The 
extra parameter v 3 introduces trigonal warping in a simi- 
lar manner as in bilayer graphene 28] , and the V4 param- 
eter produces a weak electron-hole asymmetry by adding 
the band ener gy 2 vv4p 2 /^i in both the conduction and 
valence bands |62j |. 

The Landau level spectrum in a uniform and per- 
pendicular magnetic field may be found using the Lan- 
dau gauge A = (0, Bx, 0). Then, at valley K +1 the 
operators 7r and 7rt coincide with raising and lower- 
ing operators f63l in the basis of Landau functions 



4> n {x — PyXg), such that inp n 



i{h/\ B )y/2(n + l)ip n +\, ^ ip n = -i(h/\ B )V2nip n -i, 
and 7rT ip = 0. Here A^ = yj H/(eB) is the mag- 
netic length. At valley K^, the effect of the opera- 
tors becomes n< ip n = — HH/Xb) \/2(n~+Tjip n +i, iripn — 
i(h/ As)v / 2nV'n-i > an£ l Wto — 0. In the absence of 1)3 and 
i>4, the Landau level spectrum for the Hamiltonian 1-L m m 
at the valley £ is given by 
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where L# = y/2hv 2 eB = \/2Hv/\b and we consider 
{|e|, y^Ts} -C |7i|- When a / /I, each of two lowest 
levels at n = — 1, split in valleys due to next-layer cou- 
pling 72, moving to energy corresponding to either the 
bottom of the zero-field conduction band or the top of 
the valence band. Other levels are valley degenerate in 
this approximation. 

Even the degeneracy of the n = and n = — 1 levels 
is lifted in the higher order of rs/(A7i). By applying a 



perturbation to the original 4x4 Hamiltonian, we find 
the correction to be 



Se 



n=-l ~ 
:•; — 
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(23) 



so that the splitting is proportional to B. 

The monolayerlike block for m = 0, Eq. (|TT|) . is charac- 
terized by only one parameter with a — —(N—l) / (N+l). 
The energy dispersion is given by 



4%) 



A' + a(7 2 + 75) 




A' - a (72 



75 J 



+ v 2 p 2 , (24) 



which generally has an energy gap of the width |A' — 
a(j2 — 75)| at Dirac point, and an overall energy shift of 
[A' + a (72 + 7s)]/2. In a magnetic field, Landau levels 
become 1641 



,(o) 

n>l,± 



£ (0) 



A' + a( 72 + 75) 




(25) 



Similarly to the bilayer-like block, the lowest levels at 
n = of K + and K_ split to the energy of the bottom of 
the zero-field conduction band or the top of the valence 
band, while the other levels are valley degenerate. 

The Landau levels of the different blocks are hybridized 
by the off-diagonal matrix W mm '. At the valley K + 
the wavefunction of the Landau level with index n can 
be written in form (cii^-i, C2ip n , Caipm Ciipn+i) f° r the 
bilayer-like band, and (c\ip n -\, C2ip n ) for the monolayer- 
like band. Since Wmm' is diagonal and does not include 
7r or 71^ , it only couples Landau levels of different blocks if 
they have the same index n. Furthermore, parameter 73 
mixes levels n and n + 3 within the same block, and thus 
it, together with W mm ', leads to hybridization among 
the levels of different blocks whose indices are equal in 
modulo 3. Such coupling leads to anticrossing at the in- 
tersecting point of the corresponding Landau levels. At 
the valley K_ , the wavefunction of the Landau level with 
index n becomes (ci^n+i, C2VV1, c^tj) n , c^ n -\) f° r the 
bilayer-like band, and (ciip n , c^tpn—i) f° r the monolayer- 
like band, where the index at the same position differs 
between monolayer and bilayer. As a result, the above 
rule for K + changes only for the coupling between mono- 
layer and bilayer levels, where the n-th monolayer level 
couples with the n'-th bilayer level when n—l and n' are 
equal in modulo 3. 



III. TRILAYER GRAPHENE 

According to the decomposition described previ- 
ously, the Hamiltonian in basis |</> A '° dd ' 1 ), |</ , g B '° dd ' ) ), 

|#'° dd) >, l4 B '° dd) ), l^'^, l4 B ' CV ° n) ) may be writ- 
ten in block diagonal form [42j as 

where 
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% = «(A a ) + W(l/2,0), (28) 

where A2 = v2. The off-diagonal blocks in Eq. (j2"6")l con- 
necting the monolayerlike block "Ho and the bilayerlike 
block Ti.2 are identically zero, because the basis states for 
the monolayerlike and bilayerlike blocks have different 
parity with respect to mirror reflection as argued above. 

In absence of ^3 and i>4, the low-energy Landau level 
spectrum for the bilayer-like band (m = 2) is given by 
Eq. {22} with (A,a,/3) = (V^, 1/2,0), and that for the 
monolayer-like band (m = 0) by Eq. (|25p with a = —1/2. 

The lowest Landau levels of each block e^°2 0l e i=_i an d 

(2) 

£n=(b which are degenerate in the absence of next-layer 
coupling [2!|, split according to 

a» = ^(-|) + ^(A'-|), <*> 

4£-i = i 2 2o = (i + l)|- (30) 

(2) (2) 

The two lowest levels of the bilayerlike block e n l_ 1 , e„=o, 
are weakly split by extra band parameters in accordance 
with Eq. (|2U)) . As a result, the twelve zero-energy lev- 
els split into six different energies, each of them having 
twofold spin degeneracy. 

We numerically calculate the Landau level spectrum by 
diagonalizing the original Hamiltonian Eq. ([T]) including 
all parameters. We adopt the parameter values [32[ 72 = 
-0.02 eV, 75 = 0.04 eV, A' = 0.05 eV, 7o = 3 eV, 71 = 
0.4 eV, 73 = 0.3 eV, 74 = 0.04 eV. Since the dimension of 
the Hamiltonian matrix becomes infinite in the presence 
of trigonal warping, we introduce a cut-off in the Landau 
level index, n = 100, which is enough high to obtain the 
proper low-energy spectrum 62] . 

Fig.^a) and (b) show the zero-field dispersion and the 
Landau levels plotted against the magnetic field, respec- 
tively, of ABA trilayer graphene with all the band param- 
eters included. The symbols 'M' and 'B' represent the 
monolayer block (to = 0) and the bilayer block (to = 2), 
respectively. The spectrum is composed of monolayer- 
like and bilayerlike Landau levels that are shifted rela- 
tively to each other in energy, as qualitatively described 
by the Hamiltonian decomposition above, Eqs. (|26H28I) . 
Zero energy levels are close to those obtained analytically, 
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FIG. 2: (a) Low-energy band structure and (b) Landau levels 
as a function of magnetic field of ABA-trilayer graphene. (c) 
Two-dimensional density plot of the density of states in the 
space of magnetic field and carrier density. In (b) and (c), 
numbers represent the quantized Hall conductivity in units of 
e 2 /h. 



Eqs. (|29|30|) . corresponding to the bottom of the zero- 
field conduction band or the top of the valence band. We 

(2) (2) 

also observe the weak splitting of e„Z_ 1; e„Z [indicated 
as (B,— 1), (B,0), respectively] argued above. 

Fig- El^c) shows a two-dimensional plot of the density of 
states in the space of magnetic field and the carrier den- 
sity. For simplicity, we assume that each Landau level is 
broadened into a Gaussian shape with width AE = CT b , 
[24I l65| where the constant C is taken to be 0.03. The 
dark and bright colors represent low and high density 
of states. Apart from the usual Landau fan, we observe 
a series of bright lines corresponding to the crossings of 
monolayerlike and bilayerlike Landau levels [Hf|. The 
numbers assigned to the dark regions indicate the quan- 
tized Hall conductivity in units of e 2 //i, and correspond 



6 



to those assigned to spaces between Landau levels in Fig. 
[DJb). The quantized Hall conductivity jumps by 2 only 

at the zero levels e^°2 , 4f=-i' 4i=o which split in valleys, 
while it changes in units of 4 otherwise since other levels 
are almost valley degenerate. As a result, the quantized 
Hall conductivity becomes AM (M; integer) only in the 
region between e^ 2 levels of two valleys, corresponding 
to the energy gap of the monolayer band, and also in the 



Fourlayer 



narrow regions between e^ 2 2. 
and K _ . Otherwise it takes a series of AM 



1 and e„l for each of K + 



IV. FOURLAYER GRAPHENE 

The Hamiltonian of ABA-stacked fourlayer Hamilto- 
nian, Eq. (JXJ) for N = 4, may be partially decomposed 
into two bilayerlike blocks with subsystem indices m = 1 
and m = 3: 



where 



Hi H13 
H31 H3 



Hi=H(\i) + W(-p, -p), 
H 3 = H(\ 3 ) + W(p,p), 
His = H31 = W(p/2, -p/2). 



(31) 



(32) 
(33) 
(34) 



and Ai = (-1 + VE)/2, A 3 = (1 + a/5) A and p = l/VE. 

By neglecting the inter-block mixing as well as V3 and 
t>4, the low-energy Landau level spectrum is given by 
Eq.fl22J|, with (A,a,/3) = (Ai,-p,-p) for m = 1, and 
(\3,p,p) for m — 3. The zero-energy Landau levels be- 
come 



_ ,(i) _ _72_ 

-„=_! - e„=o ^ 

c (3) . (3) = . 72_ 

-n=-X £ n=0 



(35) 
(36) 



which are valley degenerate as they should be, due to 
spatial inversion symmetry. The two lowest levels of each 

bilayer-like block, e^_ x ,e£j Q , are split due to extra pa- 
rameters as described by Eq. (|23|) . so that the sixteen 
zero-energy levels split into four energies, each of them 
with fourfold spin and valley degeneracy retained. 

Fig. [3] (a) and (b) show the zero-field band structure, 
and the Landau levels of fourlayer graphene, respectively 
computed numerically including all parameters. In (a), 
the dotted curve indicates the energy bands with the 
off-diagonal block neglected. The symbols 'b' and 'B' 
represent the light-mass bilayer (m = 1) and the heavy- 
mass bilayer block (m — 3), respectively. The Landau 
level spectrum is basically regarded as composition of 
two series of bilayerlike levels, while there are anticross- 
ings caused by H13 between the levels having the same 
index in modulo 3. The valley degeneracy is never bro- 
ken even in the full parameter model as it is protected 
by the spatial inversion symmetry of the lattice. 
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FIG. 3: Plots similar to Fig. [2] for ABA fourlayer graphene. 
In (a), the dotted curve indicates energy bands calculated by 
neglecting the off-diagonal block. 



Fig. G3c) illustrates a two-dimensional map of the den- 
sity of states calculated from the Landau levels in Fig. 
[3fb). The quantum Hall integer is always AM (M: inte- 
ger) because all the levels are valley and spin degenerate. 
Similarly to the trilayer characteristic bright lines appear 
at crossing points of Landau levels belonging to the dif- 
ferent bilayer blocks. 



V. FIVELAYER GRAPHENE 

The Hamiltonian of ABA-stacked fivelayer Hamilto- 
nian, Eq. (JXJ) for N = 5, may be partially decomposed 
into a monolayerlike block m = and two bilayerlike 
blocks with subsystem indices m = 2 and m = A: 



Hn=5 



H a Hoi^ 
H 2 
,H 40 Ha 



(37) 
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where 



Fivelayer 



(b,2+) 







2 


( 72 


) 


\VTT 


A' 




1° 


75/ 



•Ho 

H 2 =H(\ 2 )+W (0,-1/2), 
H A =H(\i) + W(2/3,l/2), 
Hoa = h1 = [W(V2/6,0)] 2> 



(38) 

(39) 
(40) 
(41) 



where Aa = 1, A4 = ^/?>, and [W]2x4 represents the upper 
half (i.e., the first two rows) of the matrix W. The block 
m = 2 is never mixed with m = and 4, due to the 
parity difference as argued previously. 

By neglecting the inter-block mixing, V3 and V4, the 
zero-energy Landau levels are given by 



£ (0) 


= '; f ( 


e n=0 




£ (2) 


- e (2) - 


fc n=-l 


— fc n=0 — 


£ (4) 


- e (4) - 


t n=-l 


— fc n=0 — 



(A' 



75 

2 



-072 



3 



72- 



(42) 
(43) 
(44) 



Similarly to the trilayer and fourlayer, the two lowest 
levels of each bilayer-like block, e^=_i, e„=o> split due to 
extra parameters in accordance with Eq. (|23l) . resulting 
in ten different zero-energy levels with spin degeneracy. 

Fig. H] (a) and (b) show the zero-field band structure, 
and the Landau levels of fivelayer graphene, respectively 
numerically computed with all the parameters. Here 
the symbols 'M', 'b' and 'B' represent the monolayer 
(m = 0), the light-mass bilayer (m = 2) and the heavy- 
mass bilayer block (m = 4), respectively. Fig. @Ic) il- 
lustrates a two-dimensional map of the density of states 
and the quantum Hall conductivity. Similarly to the tri- 
layer, the quantized Hall conductivity jumps by 2 only at 
the valley-split zero energy levels of monolayer-like and 
bilayer-like subbands, while it changes in units of 4 at 
other levels which are almost valley degenerate. As a re- 
sult, the quantized Hall conductivity becomes 4M only in 
the region between monolayer zero-levels of two valleys, 
and in the narrow gaps between —1st and 0th levels of 
each bilayer band. Otherwise it takes a series of AM + 2. 



VI. CONCLUSION 

Here, we used the effective mass approximation to ana- 
lyze the Landau level spectra of ABA-stacked multilayer 
graphenes. In general, the next-nearest layer couplings 
significantly influence the low-energy spectrum and the 
quantum Hall effect by causing energy shifts, level anti- 
crossings, and valley splitting of the low-lying Landau 
levels, as described in detail here for trilayer, fourlayer 
and fivelayer graphene. We considered a single-particle 
picture in order to provide a simple description of a broad 
range of features. Depending on sample quality, electron- 
electron interactions will also contribute to symmetry 
breaking and Landau level splitting (67T479T ]. as will in- 
terlayer asymmetry due to the presence of an external 




4 6 
Magnetic field (Tesla) 



FIG. 4: Plots similar to Fig. [2] for ABA fivelayer graphene. 
In (a), the dotted curve indicates energy bands calculated by 
neglecting the off-diagonal block. 



gate or doping |28|, |3J, |42j, |45|, |80|-|82j. Nevertheless, 
an experimental observation of features related to next- 
nearest layer couplings should be possible in high mobil- 
ity samples including suspended graphene 69l |83|-|86| or 
graphene on a boron nitride substrate (6C. l72l l73l. \8T\. 
In fact, Landau level crossings in trilayer graphene were 
recently observed [66|], and they allowed the determina- 
tion of Slonczewski-Weiss-McClure parameter values in 
remarkably close agreement with those of bulk graphite 
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